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Based  on  the  solution  of  a  simplified  canonical  problem,  the  possible 
margin  of  error  on  the  predicted  coupling  length  according  to  the  coupled 
mode  theory  as  applied  to  fiber  and  integrated  optical  guides  is  inferred. 

It  was  found  that  coupling  length  obtained  according  to  the  coupled  mode 
theory  is  usually  accurate  to  within  20%  of  the  actual  value  provided  that  the 
frequency  of  operation  is  above  the  cutoff  frequency  of  the  antisymmetric  mode 
of  the  coupled  structure. 
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Coupled  node  theory  has  been  used  extensively  in  recent  years  In  the 

1-3 

analysis  of  many  Important  problems  in  Integrated  optics  and  fiber  optics. 


It  is  therefore  of  extreme  Interest  to  learn  the  region  of  validity  for 


which  the  coupled  mode  theory  will  yield  accurate  results.  Since  exact  solu¬ 
tions  for  most  practical  problems  are  not  available,  we  shall  therefore  deal 
specifically  with  a  canonical  problem  whose  exact  solution,  as  well  as  approxi¬ 
mate  solution  based  on  the  coupled  mode  theory,  is  obtainable.  Comparison  of 
the  results  may  then  provide  an  indication  of  the  accuracy  of  the  coupled  mode 
theory  approach.  It  is  shown  that  when  the  coupled  mode  theory  is  used 
properly,  the  prediction  of  coupling  distances  is  surprisingly  accurate  even 
when  the  separation  between  the  guides  is  relatively  small.  However,  as  far 
as  the  field  configurations  of  coupled  guides  are  concerned,  significant 
Inaccuracy  is  observed  when  the  separation  distance  becomes  small. 

The  geometry  of  the  problem  is  shown  in  Fig.  1.  Two  identical  dielectric 
slab  waveguides  are  located  in  the  regions  d/2  £  |x|  £  a+d/2.  The  permittivity 
of  the  slab  guides  is  e ^  and  that  of  the  regions  outside  the  guides  is 
furthermore,  >  C2*  assumed  to  be  y^  everywhere.  A 

perfectly  conducting  plane  rovers  the  surface  z  -  0  except  for  the  region 
d/2  £  x  £  a+d/2.  In  that  aperture,  the  tangential  electric  field  is 
All  electromagnetic  field  quantities  are  Independent  of  the  coordinate  y  and 
the  time  dependence  expCjut)  is  assumed  and  suppressed. 


The  electromagnetic  field  has  components  E  ,  E  and  H  ,  which 
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where  k  *  *  up  e. ,  i  •  1,2.  Now  it  is  necessary  to  specify  that  all  field 
i  o  i 

components  have  the  same  2-dependence,  i.e.,  exp(-jfi2).  Thus  if  E^(x,z)  * 
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Ex(x)exp(-jgz)  and  similarly  for  E^  and  Hy,  we  have 
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Me  are  concerned  with  the  excitation  by  the  aperture  field  of  the  guided 
nodes  on  this  structure;  the  guided-wave  field  will  be  the  dominant  portion 
of  the  total  field  for  large  z.  Solutions  of  Eq.  (3c)  corresponding  to  the 
guided-wave  field  are 


-  •  <x  <  -a-d/z:  H  ^(x)  *  A„  explv  (x+a+d/2)]  (4a) 

—  ■*  yn  n  n 

-a  -  d/2  -  d/2:  H  ^(x)  -  B  sin  u  (x+d/2)  +  C  cos  u  (*W/2)  (4b) 

yn  n  n  n  n 

-  d/2  £  x  £  d/2:  H  ^(x)  •  D  sinh  v  x  +  E  cosh  v  x  (4c) 

■  yn  n  nn  .n 

d/2  £  x  £  a  +  d/2:  H  <4)(x)  -  F  sin  u  (x-d/2)  +  G„  cos  u  (x-d/2)  (4d) 

**  yn  n  n  o  n 

•  +  d/2  5  x  £  •:  H  *5*(x)  •  H  expl-v  (x-a-d/2)]  (4e) 

yn  n  n 


in  which  u  2  »  k  2  -  0  2*  v  2  •  0  2  -  k  2,  and  the  coefficients  A  —  H  are 
n  x  n  n  n  *  n  n 

to  be  determined.  The  *V'  subscripts  refer  to  the  nth  guided-wave  mode. 
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The  characteristic  equation  from  which  the  values  of  nay  be  deternined 

is  found  by  forcing  the  tangentical  field  components  E  and  H  to  be 

z  y 

continuous  at  the  four  dielectric  interfaces.  By  virtue  of  the  symmetry  of 

the  configuration  about  x  m  0,  the  electromagnetic  field  may  be  separated  into 

two  parts  whose  axial  electric  field  E  has  either  even  or  odd  symmetry 

z 

about  x  ■  0.  For  the  even  modes,  H  ■-  A  ,  G  «-C  ,  F  ■  B  ,  and  E  «  0; 

n  n  n  n  n  n  n 

for  the  odd  modes,  H  *  A  ,  G  *  C  ,  F  B  ,  and  D  *  0.  Thus  considering 

n  n  n  n  n  n  n 

the  two  cases  separately,  one  matches  'the  boundary  conditions  at  x  *  d/2 
and  x  *  a  +  d/2,  and  obtains 
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in  which  q  ■  e^/e^.  roots  t^>e  f*rst  factor  in  square  brackets  yield 

the  even  modes;  those  of  the  second  factor,  the  odd  modes. 

Ve  may  also  evaluate  the  coefficients  A  -  H  in  terms  of  any  single  one 

n  n 

of  them,  say  A^.  The  relationship  is  expressed  as  Bq  “  b^A^ , . . . ,  ■  h^A^  in 

which  b  -  h  are  given  by 
n  n 
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odd  modes: 
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This  completes  the  formal  analysis  of  the  guided -mode  electromagnetic  field  on 
the  dual  slab  waveguide  structure. 

We  now  consider  the  boundary  condition  at  *  -  0,  i.e.,  that  Ex(x,o)  -  Eq 
for  d/2  £  x  £  a  +  d/2  and  E^tx.o)  -  0  elsewhere.  Equating  the  excitation 
field  to  the  total  field  E^  for  x  jj  0  at  t  •  0  yields 

E0[lH*-4/2>  -  lKx-.-d/2)]-  £  3^0  V«>  <8> 

n 

+  (radiated  field)  | 
x*0 

where  U  is  the  unit  step  function.  The  sum  is  taken  over  the  guided  modes  for 
which  Bn  la  a  proper  root*  of  Eq.  (5).  By  virtue  of  the  orthogonality  relation 
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The  integrals  in  Eq.  (10)  are  easily  evaluated  and  yield  a  relation  between  E 

o 

and  A  .  We  find  that 
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in  which  the  lover  (-)  sign  is  taken  for  the  even  modes »  and  the  upper  (+) 

2  e2  o2 

sign  for  the  odd  nodes.  4^  is  either  4^  or  4”'  as’ appropriate,  as  are 
fQ  and  gfi.  This  completes  the  formal  solution  for  the  excitation  of  the  guided- 
wave  field  by  the  aperture  electric  field  at  t  ■  0.  Here  we  have  assumed 
that  the  contribution  of  the  exciting  field  to  the  radiated  field  is  negli¬ 
gible  as  compared  with  the  guided  field. 

We  now  consider  the  special  case  which  is  of  particular  interest  for  the 
purposes  of  this  paper,  that  in  which  only  a  single  node  is  possible  on  each 
of  the  dielectric  slabs  in  isolation.  We  shall  further  require  that  the  * 
loves t-order  even  and  odd  nodes  will  propagate  on  the  dual-guide  structure. 
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Thus, 


2a  „ 

— r  <  ua  can  ua  <  • 

qa 

where  0  <  ua  <  e/2.  We  now  evaluate  the  r-component  of  the  Poynting  vector 
in  each  of  the  slab  waveguides  under  the  assumption  that  only  these  lowest- 
order  modes  exist;  their  propagation  constants  are  denoted  6  and  6  , 

€  O 

and  the  associated  fields  H  (x)  are  denoted  H  (x)  and  H  (x)  respectively. 

y  ,  ye  yo 

In  the  "upper”  guide  (d/2  <  x  <  a+d/2)  or  "guide  A",  we  have 
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In  the  "lower"  guide,  or  "guide  B", 
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Ve  may  now  calculate  the  total  power  per  unit  width  carried  in  each  of  the 
(  + 

•lab  waveguides  by  integrating  from  x  •  d/2  to  x  ■  e+d/2  and  from 
x  a-d/2  to  x  "-d/2.  Denoting  these  powers  by  P^  and  P^,  we  obtain. 


^The  power  carried  outside  the  slab  guides  is  not  included  in  the  calculation. 
Thus  PA  and  Pg  will  not  reduce  to  the  total  power  carried  in  a  single  slab 
waveguide  mode  in  the  limit  d  •. 
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assuming  that  c^»  tj,  B#.  and  #e  arc  purely  real, 
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Now  ve  are  in  a  position  to  compute  ev»ctly  how  the  input  power  is 

transferred  from  one  guide  to  the  other.  (The  only  assumption  is  that  no 

radiated  power  is  present.)  We  must  first  determine  the  propagation 

constants  of  the  guided  modes  along  the  slabs  from  Eq.  (5).  Results  for 

the  two  lowest  order  modes  are  displayed  in  Fig.  2  for  the  symmetric 

modes  and  in  Fig.  3  for  the  anti-aynsetric  modes.  In  these  figures,  the 

normalized  propagation  constant  B/ko  is  plotted  against  the  normalized 

thickness  of  the  slab  k  a  for  various  values  of  the  normalized  separation 

o 

oi  the  slabs  k  d.  The  constant  k  ,  the  free-space  wave  number,  is  defined 
o  o 

as  «tVp  c  .  It  can  be  seen  that  there  exists  no  cutoff  frequency  for  the 
o  o 

lowest  order  symmetric  mode  while  a  cutoff  frequency  does  exist  for  the  lowest 
order  anti-sysnetric  mode.  •  Of  course,  all  higher  order  symnetric  and  anti¬ 
symmetric  modes  possess  cutoff  frequencies. 

To  illustrate  how  the  guided  power  exchanges  between  the  two  guiding 
structures.  Figs.  4-5  are  Introduced.  The  operating  frequency  is  so  selected 
that  only  the  dominant  symmetric  mode  and  the  lowest  order  antisymmetric  mode 
may  exist  along  the  guiding  structure.  It  is  interesting  to  note  that  although 
the  initial  exciting  field  exists  only  at  the  entrance  of  the  core  region  of 
guide  A,  according  to  our  computed  results  there  exists  a  small  amount  of  guided 
power  in  guide  B.  The  reason  for  this  is  that  to  satisfy  the  initial  given 
field  configuration  at  the  entrance  of  the  guiding  structure,  radiation  mode 
as  well  as  guided  modes  must  be  taken  into  account.  Since  we  have  a  priori 
ignored  the  radiation  mode  in  our  calculation,  ve  can  only  satisfy  approximately 


the  given  initial  field.  The  quantity  it  *  ■  0  represents  the  power  calcu¬ 
lated  from  the  field  that  extends  from  guide  A  into  the  core  region  of  guide  B. 
It  can  also  be  seen  from  these  figures  that  guided  power  exchanges  from  one 
guide  to  the  other  in  a  periodic  fashion  as  expected.  The  distance  for  which 
maximum  guided  power  is  transferred  from  one  guide  to  the  other  is  called  the 
coupling  length.  It  is  noted  that  the  coupling  length  becomes  shorter  as  the 

separation  distance  between  the  guides  is  shortened.  To  correlate  the  maximum 

* 

power  contained  in  the  core  region  of  guide  A  with  the  normalized  separation 

distance  k^d,  ve  have  performed  the  computation  at  the  entrance  of  the  coupled 

guide.  Results  are  shown  in  Fig.  6.  It  can  be  seen  that  the  normalized  maximum 

power  in  the  core  region  of  guide  A  varies  in  a  rather  unusual  fashion  for  small 

separation  distances.  At  large  k^d,  the  maximum  power  contained  in  the  core 

region  approaches  that  for  the  case  of  an  Isolated  slab  guide,  as  expected. 

To  further  understand  the  behavior  of  the  transverse  fields  in  the  guides  when 

the  separation  distance  is  small,  we  have  plotted  in  Fig.  7  the  quantity  |E^ | 

vs.  the  transverse  distance.  The  complexity  of  the  evolution  of  the  transverse 

electric  field  as  the  separation  distances  are  changed  indicates  the  complex 

nature  of  the  curves  in  Fig.  6  when  the  separation  distances  are  small. 

Recall  that  the  primary  purpose  of  our  investigation  is  to  determine  how 

accurate  the  coupled  node  theory  is  in  its  treatment  of  the  coupled  dielectric 

waveguide  problem.  Ve  have  carried  out  the  cases  treated  above  according  to 

1  2 

the  coupled  mode  theory  described  by  Marcuse  and  by  MELC  researchers.  The 
"exact"  normal  mode  results  are  then  compared  with  those  obtained  according 
to  the  various  coupled -mode  theories.  Displayed  in  Fig.  8  are  the  curves  for 
tits  coupling  length  as  a  function  of  the  normalized  separation  distance  be¬ 
tween  two  parallel  dielectric  guides  as  shown  for  two  different  k  a  values  in 
Fig.  1.  One  notes  that  as  the  separation  distance  is  increased,  the  agree- 


sent  between  tbe  results  based  on  coupled  mode  theories  and  our  normal  mode 

results  becomes  better,  end  that  closer  agreement  Is  obtained  for  larger 

fcoa  values  or  when  more  power  is  confined  within  the  core  of  the  guide.  This 

is  because  large  koa  corresponds  to  more  tightly-bounded  field;  hence,  the 

coupling  field  is  weaker  and  the  coupled  mode  theory  which  is  based  on 

perturbation  concept  tends  to  be  SK>re  accurate.  It  is  worthwhile  to  point 

out  that  for  certain  combinations  of  k  a  and  k  d  values  such  as  for  k  a  *  1 

o  o  o 

and  k^d  £  7.0  only  the  dominant  syuetrlc  mode  exists,  so  according  to  the 
exact  normal  mode  theory,  no  back  and  forth  exchange  of  propagating  power  takes 
place  between  guide  A  and  guide  B.  On  the  other  hand  the  approximate  coupled- 
mode  theory  continues  to  predict  the  power  exchange  phenomenon.  Another  way 
of  expressing  the  differences  for  the  results  based  on  different  theories  is 
shown  in  Tig.  9  where  the  percent  differences  between  different  coupled-mode 
theories  and  the  normal  mode  theory  are  plotted  against  the  normalized 
separation  distances. 

It  can  be  seen  that  the  coupled  mode  theory  is  surprisingly  good  (within 
20Z)  in  predicting  the  coupling  length  of  two  parallel  dielectric  slab  guides 
even  when  the  separation  distance  Is  relatively  small  and  the  confinement  of 
guided  power  is  relatively  weak.  Extrapolating  our  present  results  to  other 
geometries  involving  optical  fibers  or  integrated  optical  guides,  It  is 
Inferred  that  the  coupling  distances  predicted  according  to  the  coupled  mode 
theory  are  accurate  to  within  20X  of  the  actual  values  if  .the  symmetric  and 
antisymetrlc  modes  are  both  above  cutoff.  Finally ,  it  should  be  noted  that 
when  the  aeparatlon  distance  is  small  the  transverse  field  configurations  of 
the  coupled  guides  (see  Fig.  7)  are  significantly  different  than  those  assumed 
in  the  coupled  mode  theory. 
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Geometry  of  the  canonical  problem. 

Normalized  propagation  constants  as  a  function  of 
normalized  frequencies*  The  lowest  order  mode  n  »  1 
has  zero  cutoff  frequency. 

Normalized  propagation  constants  as  a  function  of 
normalized  frequencies.  All  antisymmetric  modes  have 
cutoff  frequencies . 

Normalized  power  in  the  core  region  of  guides  A  and  B 
as  a  function  of  the  normalized  longitudinal  distance. 

The  coupling  length  (defined  as  the  length  for  which 
complete  exchange  of  power  in  the  cores  of  Guide  A  and 
Guide  B  occurs)  is  longer  for  less  tightly  bounded-fields. 

Normalized  power  in  the  core  region  of  guides  A  and  B  as 
a  function  of  the  normalized  longitudinal  distance.  P 
is  the  power  in  the  core  region  of  guide  A  or  guide  B  as 
appropriate.  P^,  is  the  total  guided  power. 

Maximum  normalized  power  in  the  core  region  of  guide  A  as 
a  function  of  separation  distance  of  the  two  guides. 

Transverse  electric  field  distribution  across  the  two 
coupled  guides. 

Normalized  coupling  length  as  a  function  of  the  normalized 

separation  distance.  Note  that  the  coupling  length  ceases 

to  exist  for  k  d  <7.0  when  k  a  «  1.0  according  to  the  normal 
,  -  o  o 

node  theory. 

Percent  coupling  length  differences  between  normal  mode 
theory  and  coupled  mode  theory  as  a  function  of  normalized 
separation  distance. 
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